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Abstract 

We study a classical reparametrization-invariant system, in which "time" is not a priori 
defined. It consists of a nonrelativistic particle moving in five dimensions, two of which are 
compactified to form a torus. There, assuming a suitable potential, the internal motion 
is ergodic or more strongly irregular. We consider quasi-local observables which measure 
the system's "change" in a coarse-grained way. Based on this, we construct a statistical 
timelike parameter, particularly with the help of maximum entropy method and Fisher- 
Rao information metric. The emergent reparametrization-invariant "time" does not run 
smoothly but is simply related to the proper time on the average. For sufficiently low 
energy, the external motion is then described by a unitary quantum mechanical evolution 
in accordance with the Schrodinger equation. 

PACS: 04.20.-q, 04.60.Kz, 03.65.Ta, 05.20.-y 

1 Introduction 



Motivated by attempts to quantize gravity, based on the classical theory of General Relativity, 
there has recently been interest in the quantization of "timeless" reparametrization-invariant 
systems, see, for example, Refs. |L], 0, [| |J and further references therein. 

Presently, we begin with the study of a classical system and address the question, whether 
local observables can be found which allow to characterize the evolution in a gauge invariant 
way. In previous work it has always been assumed that the global features of the trajectories 
are accessible to the observer, which makes it possible, in principle, to express the evolution 
of an arbitrarily selected degree of freedom "relationally" in terms of others || Thereby 
the Hamiltonian and possibly additional constraints have been eliminated in favour of Rovelli's 
"evolving constants of motion" 

In distinction, we presently attempt to characterize the evolution by invariant quasi-local 
statistical properties of the ergodic internal "clock" motion. Heuristically, we assume "time is 
change" and try to quantify the former in terms of measurements of the latter. Our results 
indicate that a 'deparametrized' time evolution can be constructed based on coarse-graining 
localized observations in a classical reparametrization-invariant system which is ergodic. 
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We remark that certain forms of globally incomplete statistical knowledge about a classical 
system lead to its effective quantization locally J7|. This points towards a deterministic origin 
of quantization and certainly raises further interesting questions about the relation to the 
problem at hand. Indeed we shall find that the external motion of the particle is described by 
a (discrete-time) Schrodinger evolution. 

Let us consider a fivedimensional model of a "timeless" nonrelativistic particle with the 
action: 

S= [dtL , (1) 



where the Lagrangian is defined by: 



L = 



2A 



{d^-r^d^f-r 2 ^) 



A 
2 



2r>2; 



W2 



-E 



(2) 



Here A stands for an arbitrary "lapse" function of the parameter t, q e R 3 denotes an ordinary 
vector, and 0i and 02 denote the angular variables corresponding to the toroidally compactified 
dimensions with radius r, respectively; u> 2 , Q 2 are angular velocity squared coupling parame- 
ters. The parameter E fixes the energy difference between the "external" and the compactified 
"internal" degrees of freedom. 

Suitably redefining the units of length and energy, we set r = 1 henceforth. The notation 
indicates that the corresponding terms in L are periodically continued: 



n 



E [n, n + 1[ 



(3) 



for any integer n. Thus, a ratchet type potential results in the 01,2-plane. Alternatively, we 
may consider the angular variables to be normalized to the square [0, l[x[0, 1[, of which the 
opposite boundaries are identified, thus describing the surface of a torus with main radii 1/27T. 

The relative minus sign between the kinetic terms allows the energy of the external motion 
to be unbounded from above. A similar feature is common to simplified cosmological models 
in accordance with General Relativity, as discussed, for example, in Refs. || £| and references 
therein. Furthermore, we choose the potential in the compactified dimensions to be unstable, 
in order to generate a chaotic internal motion which is ergodic and possibly more strongly 
irregular. 

Varying the action with respect to q, 0i, and 02, we obtain the equations of motion: 



\d t (\d t q) 
1 1 

-rd t {-d t (j) 2 ) 



, 

(^ 2 1 + fi 2 2 )(i-E^i-^) ) 

n 



(4) 
(5) 

(6) 



where the singular terms arise due to the discontinuities of the potential. Variation of S with 
respect to A yields the constraint: 



1 

2P 



(d t q) 2 - {dt^f - {d t <h 



+ 2 



W 2 (0 2 + 2 2 ) + 2fi 2 0!0 2 



E 



, 



(7) 



which will be recognized as a constraint on the Hamiltonian momentarily. 
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The canonical momenta are defined as usual: 

dL 1 

P = 



d{d t q) A 
dL I 



7Tl 



7T 2 = 



didtfa) A 
<9L 1 



d{d t (h) A 
Incorporating these, we obtain the Hamiltonian: 

H = p ■ d t q + 7Ti9 t 0i + 7ri<9 t 2 - £ 

A 
2 



d t <p2 ■ 



V 



AC . 



In terms of C, Eq. (0) implies the primary constraint: 

C = , 

which is a weak equality in the sense of Dirac's formalism of constraint systems 
quently, the Hamiltonian presents a weak constraint, H = 0. 

Finally, in Hamiltonian form, the equations of motion read: 

d t p = , 

n 

d t ir 2 = -\(u% + n%)(l-J2S(<f>2-n) ) . 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 
Conse- 



(14) 
(15) 

(16) 



Employing these equations and the definition of the constraint C in Eq. (0), it follows by 
explicit calculation that it does not evolve. Consistently, this is also obtained by: 



d t C = {C, H} = \~ {H, H} = 



(17) 



employing the Poisson bracket notation, {A, B} = J2(dQAdpB — dpAdqB), where a sum over 
all coordinates and canonical momenta, represented by Q and P, is understood. Therefore, no 
secondary constraints exist in our model. 

We conclude that the system has four physical degrees of freedom. Its extended phase 
space, cf. Ref. [0, is tendimensional, corresponding to the Lagrangian variables in Eq. @ and 
the associated canonical momenta. It is, however, reduced to a ninedimensional surface by 
the constraint. Since the physical phase space is eightdimensional, there must be one "gauge" 
degree of freedom, which is related to the reparametrization invariance. The study of this 
gauge symmetry and its consequences will be performed in Section 2 , while the gauge invariant 
description of the evolution will be developed in Section 3 . In Section 4 we demonstrate that 
under certain conditions the external motion of the particle can be mapped onto an evolution 
according to the Schrodinger equation. We conclude with a brief discussion. 
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2 Gauge invariance and observables 

We observe indeed that the action is invariant under the set of gauge transformations: 

f = /(f) , x(t) = x'(t') , A(f)^ = A'(0 , (18) 

where x G {q, 0i,02}- Corresponding infinitesimal transformations are generated by: 

8t = t-t' = e(t') , (19) 
where e is infinitesimal. This yields immediately: 

5x = x{t') - x'(t') = -e{t')d t ,x{t') , (20) 

5X = A(0-A / (0 = -5f[e(* / )A(0] • ( 21 ) 
Employing the definitions of the canonical momenta, Eqs. fl8|)— fllOD, we obtain from Eq. (|20|) : 

5q = -e\p , (22) 

50i,2 = eA7Ti i2 . (23) 



Similarly, with the help of Eqs. (|14[)-([16|), we obtain: 

6p = , (24) 
S7n = eX(u% + Q 2 ^ 2 )(l-J2^<Pi-n) ) , (25) 

n 

5n 2 = eX(to% + n%)(l-Y,S((p2-n)) . (26) 

n 

Comparing with the equations of motion, we see that the evolution of the coordinates and 
momenta is generated by the gauge transformations. 

Obviously, the three-momentum p is conserved, in accordance with translation invariance. 
Its components present three gauge invariant observables which may serve as coordinates of 
the physical phase space. Another invariant is provided by the internal contribution to the 
constraint: 

C mt = p 2 - E - 2C , (27) 

cf. Eqs. (|i"2"D-([T3|); since the constraint does not evolve, i.e. is invariant, as well as p 2 , this also 
holds for Ci nt . However, being related via the constraint, Ci nt is not independent of p. The 
angular momentum suggests itself as a further observable, 

J = qxp, (28) 

due to rotational invariance. Clearly, 5J = 5qxp + qx5p = 0. 

As expected, the external motion plays a rather passive role in our system, since it can be 
almost completely described in terms of the conserved linear and angular momenta. It only 
contributes with its kinetic energy to the constraint C . However, note that p and J together 
determine only the two constant components of q transverse to p, i.e. q±. Therefore, in order 
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to predict the evolution of the coordinate q, we need to construct a reference "clock" , such that 
the longitudinal component of Eq. (|j) is well determined, when integrated: 

q\\(t) = 4 0)+ pJ» dt ' X W ■ ( 29 ) 

This is one of the objectives of the following section. 

Furthermore, we remark that we so far obtained five gauge invariant coordinates (p, q±) 
for the eight dimensional physical phase space plus one constraint (C, or C^), which allows 
to eliminate one more of the remaining four internal variables (4>i t 2, ^1,2) ■ Due to the intrinsic 
nonlinearity of the gauge tranformations, Eqs. (|22|)-(p6|), we are unable to find further invariants 
in the general, interacting case. This motivates our attempt to construct other, statistical 
observables, in order to obtain an invariant, even if coarse-grained description of the physical 
phase space. 



3 Ergodicity and coarse-grained "timing" 

In order to proceed, we make the crucial assumption that our model forms an ergodic system ||. 
We also restrict the allowed lapse functions A to be (strictly) positive, thus avoiding trajectories 
which trace themselves backwards (or stall) 0. 

In the following, we will consider quantities -AT[0i )2 ] (cj) which are functionals only of the 
trajectory determined by the coordinates (f>i^ and possibly depend parametrically on further 
constraints involving only them. While some explicit examples will be studied in detail, gen- 
erally speaking, we have quantities in mind which reflect properties of Poincare sections of 
the full trajectories. Naturally, we require the constraints q to transform as the coordinates 
under the gauge transformations QT%D, i.e. c' { = C{. It follows that such quantities are gauge 
invariant, since they depend only on geometrical properties of a path and related constraints. 
Being independent of the momenta, they do not depend on how the trajectory is parametrized: 
iV[0' 12 ](C') = N[<f>i t 2](ci) . Thus, they qualify as coarse-grained observables characterizing the 
internal motion. 

Our aim is to construct a timelike variable based on such observables. In the following first 
subsection we do this based on the idea that the geometric path length covered by the system 
evolving from an initial to a final state is an invariant measure of the "time" that passed. The 
crucial point is that this measure can be inferred in an ergodic system approximately from 
coarse-graining localized observations, provided we understand the dynamics sufficiently well. 

In the second subsection, however, considering the interacting nonlinear system, we gen- 
eralize this approach, in order to extract a "time" from quasi-local measures of the "change" 
occurring while the system evolves. In particular, we will employ a maximum-entropy method 
together with the Fisher-Rao information metric, in order to characterize the distance, i.e. the 
"time" passed, between evolving probability distributions. We point out that our approach 
is somehow orthogonal to the one of Ref. l lQfl , although we make use of its formalism. There 



the author launches the ambitious project to derive dynamics from rules of inference and the 
maximum-entropy principle in particular. We instead assume the reparametrization-invariant 
dynamics to be given and construct a pertinent notion of "time" 



6 



H.-Th. Elze and O. Schipper 



3.1 Free internal motion on the torus 

Illustrating our approach, we begin with the noninteracting case, i.e., with uo 2 = Q 2 = 0. 
Even without interactions, the internal motion is ergodic for almost all initial conditions. In 
particular, if the ratio of the two independent angular velocities is not a rational number, then, 
for sufficiently late times t, the trajectory will come arbitrarily close to any point on the surface 
of the torus. This is easily seen in its parametric representation, Eq. (|3l|) below. In Fig. 1 we 
show a typical example. 

The equations of motion immediately yield the solutions: 

Mt) = tdt'Xit') , (30) 

Jo 

where ($\ and nf\ denote the initial coordinates and momenta, respectively; of course, 7r 12 (t) = 
7r[°2. Eliminating the integral of the lapse function, we obtain, for example, 02 in terms of fa: 

(o) 

faifa)=<pf -%(0l o) -0i) • (si) 

Following Refs. [|, this is a gauge invariant "relational" description of the motion: 2 (s) 
gives the value of coordinate 02, when ('the time is such that') 0i has the value s. In this way, 
0i may serve as a time variable, even if not a unique one. 

However, similarly as in the models studied previously, there is additional information about 
the full path, which is necessary to complete this description. In the present case, the solution 
(|5ID presents a line in the 01^-plane, unwrapping the motion which multiply covers the torus. 
Thus, when folding it back onto the torus, one has to keep track of which unit square in the 
plane a respective piece is coming from. This can be labelled by two integers ni j2 , which may 
be interpreted as the winding numbers characterizing additionally a given point (0i,02) on 
the path. Obviously, this presents highly nonlocal (topological) information, which will be 
unavailable for a local observer under more realistic circumstances, such as in the presence of 
nonlinear interactions. 

Therefore, we turn to statistical measures of the motion on the torus. Consider the "incident 
number" lf(L) which counts the number of times that a given trajectory (0i(t), fait)), cuts the 
0!-axis in the neighbourhood of the observer, i.e., with fa £ [0,5 < 1] and 2 = 0, subject to 
the constraint that the total path length equals L, taking into account the wrapping around the 
torus. Due to ergodicity and assumed nonnegative lapse functions, I\ is a stepwise increasing 
function of t, while fapkt) £ [0, 1[ are of sawtooth type, with details depending on the particular 
parametrization. The reparametrization- invariant value of I\, however, only depends on the 
path length; thus, it implicitly involves nonlocal information. 

The incident number I\ presents a very simple example of the coarse-grained observables 
discussed before. There exists an unlimited number of different such counting observables. 
The more of them we introduce and measure, the more detailed will be our reparametrization- 
invariant description of the internal motion. We will also use the incident number I 2 which is 
defined like I\, however, with the roles of coordinates "1" and "2" exchanged. 

Furthermore, we may consider one incident number as a parametric function of the other: 



J 2 i(fc) = max I 2 , with ifiL) = k , 



(32) 
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where S is fixed, while the total path length L is the implicit common variable. Note that I21 is 
unique, since we take the respective maximum of I 2 , which may increase while I\ stays constant 
temporarily 

However, the constraint on a given path length is irrelevant for the actual values assigned 
to I\ and O ne operationally determines them by counting the localized incidents, with 
0i,2(^) ^ [0,5], and recording one incident number as a function of the other; no knowledge of 
the path length is required. 

Nevertheless, these reparametrization-invariant numbers can be used to determine the cor- 
responding path length by a statistical consideration, since we have sufficient knowledge of the 
dynamics of the system in the present (integrable) case. 

In the absence of interactions the path is composed of straight line segments. We denote 
the average length of these segments by (/). It can be calculated easily due to ergodicity: 

_ [(vri 0) ) 2 + (4 0) ) 2 ] 1/2 = [2E mt ]^ 



where E int denotes the conserved internal energy, cf. Eqs. fll2|) and (P7|). Here we employed 
Eq. (pip and averaged the straight-line paths with the asymptotically uniform density over the 
square [0, l[x [0, 1[, choosing coordinates such that tc^I > 0. 

Now, each increment of 1\ or /21 (i-e. I2) by one unit corresponds to the completion of one 
line segment. Due to the finite window size 5 < 1, however, only a corresponding fraction of 
all incidents happening on both coordinate axes will be recorded on the average. Correcting 
for this, invoking ergodicity as before, the path length L which leads to the measured incident 
numbers is simply obtained by: 

where ds is along the curve of the function l2i{s). Practically, what appears on the right-hand 
side is the total registered incident number (7 = I± + I 2 ) at a certain instant. 

In this way, we obtain a measure of the "time" interval T during which N incidents occurred: 

T = L/[2£ m *] 1/2 , (35) 

i.e., dividing the path length by the constant velocity. 

The point of this rather trivial example is that nowhere do we make use of the time 
parametrizing the evolution nor of the generally unknown path length. Rather, we derive 
T from reparametrization-invariant localized measurements. 

Not surprisingly, the resulting "time" will not run smoothly, due to the coarse-grained 
description of the internal motion: as if we were reading an analog "clock" under a stroboscopic 
light. This is precisely the role of a Poincare section with respect to the increasing invariant 
path length of an evolving trajectory. 

In order to illustrate the behavior of the "time" T, it is convenient to introduce the fictitious 
proper time (function): 

rt dt'X(t') . (36) 


Then, keeping the notation as simple as possible, with x(t) = x(t) for x G {q,p, 0i,2, ^1,2}, we 
obtain: dtx(t) = X(t)d T x(r). Applying this transformation to the definition of the canonical 
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momenta and the equations of motion, the lapse function can be eliminated; this replaces t 
by t and A by 1 in Eqs. (|8"l)-(|T0l) and Eqs. (|14j)-(|16|), respectively The resulting equations are 
reparametrization- invariant. Solutions of these equations are to be interpreted 'physically' by 
introducing the inverse function t{r) and x(t) = x(r(t)) 

Integrating the free internal motion with respect to the proper time, which replaces the 
integral in Eq. ( p0[) by r, we show T as a function of r in Fig. 2a) and Fig. 2b). For the two 
runs differing in the total (computer) time r, we find that after a short while, i.e. already at 
low incident numbers, the constructed physical "time" T approximates qualitatively well the 
proper time r. The fluctuations on top of the observed linear dependence naturally reflect the 
stochastic internal motion. The fact that the slopes are consistently smaller than one can be 
attributed to the bias towards longer-than-average pieces of trajectory, which is introduced by 
measuring the incidents close to the origin and extrapolating from there, see Eqs. ([33|) and (fM]). 

Finally, employing the reparametrization-invariant "time", which runs approximately par- 
allel to the unphysical proper time, T w kt (k constant), we succeed to describe the external 
motion with respect to the physical internal "clock" constructed here. Recalling that the ex- 
ternal motion is given by Eq. (|29|) , we obtain: gj|(r) = + rp ps qS ' + K~ x Tp. 

Interestingly, the coarse-grained 'jumpiness' of T introduces a corresponding stochastic com- 
ponent into the external motion. This would be recognized, if global relational data or results of 
increasingly fine-grained local measurements were available, with which to compare the change 
of q». We will further study the consequences of this stroboscopic effect in Section 4. 



3.2 Quasi-local analysis of "time is change" with interactions 

Now we consider the case of the internal motion on the torus with the interactions (oc u 2 , Q 2 ) 
turned on. The interactions additionally mix the trajectories in phase space, since the potentials 
are of inverted oscillator type. Slightly simplifying the ensuing calculations, we set u 2 = Q 2 
henceforth. This leaves the motion parallel to the diagonal 0i = 02 exponentially unstable, 
while the orthogonal motion is free. A typical trajectory for the case of still relatively weak 
interaction is shown in Fig. 3 . 

The microstates of the internal part of the system can be described by the phase space 
coordinates x = (0i, 02; 7Ti, ^2)- However, in general, we will not be able to follow the de- 
terministic evolution through sequences of microstates for any sufficiently complex nonlinear 
system. Therefore, we develop a coarse-graining statistical approach based on probability dis- 
tributions. Let P(x)dx denote the expected number of microstates in the volume element 
dx at x. Then, we may characterize macrostates of the system by giving their coordinates 
Q\ i = 1, . . . , n in the n-dimensional state space: 

Q i = (Q l ) = JdxP{x)Q\x) , (37) 
with {Q 1 } denoting the relevant set of observables. It is assumed that all information necessary 



to answer our particular questions about the system is encoded in these observables [10 



The conserved internal energy, see Eq. (|27|) , presents an important observable: 

H int (x) = 7T 2 + TX 2 - CU 2 (0x + 2 ) 2 , (31 



where we absorbed an inessential factor 2 into the definition for convenience. 
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Then, the 'prior' distribution P c (x\H int ) which best reflects our mostly lacking information 
about the state of the system, given the conserved energy, is obtained by maximizing the 
entropy: 

3[P] = - [dx P(x) lnP(x) , (39) 



subject to the constraint (Hi nt ) = H int , with this constant being fixed by the initial conditions. 
The result is the canonical distribution: 

P c (x\H int ) = Z-V^-W , (40) 

with the partition function Z c and Lagrange multiplier /3, respectively, to be calculated from: 

Z c = Jdxe-^W , H mt = -—\nZ c , (41) 

where J dx = /f^, d7i"i d^ Jq d<f>i Jq d(f>i. Thus we recover expressions which are familiar 
from statistical mechanics. 

Straightforward calculation yields the partition function: 

Z c = ^Z c , conf = ^J(l,l) , (42) 
where Z c ^ con f denotes its configurational factor. It is given here in terms of the integral: 

(5lo 2 Jo 



J(a,b) = — dsj ds 2 e {s ^ (43) 



1 7T 



a + b) erfi((a + b) J f3u 2 ) - a erfi{aJ(3u 2 ) - b erft{bJ /3u 2 ) 



2\ pu 2 

1 _ ^(a+b) 2 ^ 2 _ e a 2 ^ 2 _ e b 2 fiu> 2 + ^ ^ 

= ab (l + (3u 2 \^a 2 + ^ab+ h) 2 } 

+\{^ 2 f[\a i + \a% + \a 2 b 2 + \ab 3 + \b'} + 0[(/^ 2 ) 3 ]) , (45) 

involving the imaginary error function, erfi(x) = 1ix~ x l 2 Jq ds exp(s 2 ). We will make more use 
of this integral in the following. 

For example, incorporating the small-/? expansion, we calculate the relation between Hi nt 
and (3: 

^7T = A---— (3u; 2 + 0[((3lu 2 ) 2 ] , (46) 
u 2 (3uj 2 6 180 v ; 

using Eqs. ([11]), fl42]), and (45). Here the error in comparison with the exact result rapidly 

decreases with energy and is less than 5% for H int /uj 2 > 1. We are presently interested in the 

positive energy regime, since only there the trajectories can explore all of the torus surface, cf. 

Eq.©. 

In order to improve the prior distribution P c with the help of quasi-local measurements, we 
again consider simple observables for illustration. We define two "window functions" : 



I 1 {x) = 0(0i - e)S(5 - 0i)e(0 2 )6(e - 2 ) , (47) 
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with 0<e^C<5<l;Z 2 is denned by the analogous expression with (pi and 2 exchanged. Thus 
X\ and X 2 project out small rectangles along the <px~ and 2 -axis, respectively, which do not 
overlap. Here we explicitly introduced the small but finite window width e, which in any case 
is necessary for a correct counting of incidents in numerical simulations with a finite resolution. 
In the following we will adapt to the present case the measurement of incident numbers. 

We decompose the prior distribution with respect to the windows determined by X\ and X 2 . 

P c (x\H mt ) = [X 1 {x)+J 2 {x)\P c {x\H int ) + [l-l 1 {x)-X 2 {x)\P c {x\Hint) 

= P' w {x\H int ) + P w {x\H int ) . (48) 

While P w describes their complement, the distribution P' w is the one which is related to mea- 
surements within the windows. Its normalized counterpart P w is: 

p , ,„ v [l 1 (x)+l 2 (x)]P c (x\H mt ) (3 [J 1 (x)+X 2 (x)]e-^W 

r wK x\n mt} _ [Xi{x)+Ux)]PMHint) 2n J(5)£) _ J(CjC) > I y J 

employing the configurational integral of Eq. (|43|) . 

Then, the distribution P w (x\Hi nt ; (Xi)) which best reflects the information contained in the 
prior distribution P w (x\Hi nt ) and in the data from measurements of the window functions is 
obtained by maximizing the entropy [ |T0|j : 

S [P] ee f dx [X 1 (x)+X 2 (x)]P(x)\n ^ , (50) 

tii ^ | int ) 

subject to the constraints: 

(£> = T^T- • (51) 

-<1 + J 2 

Here we determine the average of the incident functions, i.e. the (total) probabilities of observ- 
ing an incident in the respective windows, in terms of the measured incident numbers. This 
obviously presents a crude coarse-graining. An improved description is obtained, for example, 
by binning the incident numbers with respect to the main axis of each window. Many, more 
detailed measurements can be envisaged, but the simplest ones will suffice here. 

It is straightforward to show that this procedure yields a grand-canonical distribution: 

P w (x\H int ; (X» = Z-\X x {x) + X 2 ( x )] e -^-(-)-^W-^W , (52) 
where, in this case, the partition function and Lagrange multipliers are determined by: 

dx [Ji(ar) + j 2 ( x )] e -/3W-iW-a 1 i 1 (x)- a2 x 2 (xO 

= Jdx [X 1 (x)X 1 +X 2 (x)\ 2 ]e~ m ^ x) , (53) 
introducing the fugacities Xi = exp(— a«) |i=i, 2 , and: 

(Ti) = -— In Z = Xi— In Z | i=1 , 2 , (54) 
ooti oXi 

together with the constraints (0). Here j3 is considered to be a known feature of the prior 
distribution P c , previously determined in Eq. (f46|). 
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The present situation differs in an important way from the usual one in statistical mechan- 
ics, where H int and would all correspond to conserved quantities and be treated on an 
equal footing. Presently, the (window functions related to the) incident numbers are evolving 
quantities which we measure, in order to learn about the change occurring in the system. 

The grand-canonical partition function can be calculated directly, resulting in: 



7T 



7T 



Z eonf = 7i (\i + \ i )[J(6,e)- J(e,e)] 



where Z con f is the configurational factor of this partition function. Furthermore, we obtain 

A, 



(55) 



Ai + A 2 



\i=l,2 j 



(56) 



which implies Ai/A 2 = (Zi}/(X 2 ). We set A, = C(li) (1=1,2, with a common (undetermined) 
constant C. Incorporating these results, the distribution follows: 



P w (x\Hi n ti (Zi)) 



{3 [(^Mx) + (l 2 )l 2 (x)} mi 
vr J(5,e)- J(e,e) 



(57) 



using (Xi + 1 2 ) = 1 and cancelling the common factor C. 

Finally, re- normalizing P w (x\H int ; (Z;)) and using the resulting distribution in place of P' w , 
cf. Eq. (fP|), we obtain the properly normalized distribution for the whole phase space: 

P(x\H mt ; (Z,» = 2[(J 1 )J 1 (x) + (l 2 )l 2 (x)]P c (x\H mt ) + [1 - l 1 {x) - l 2 (x)]P c (x\H int ) , (58) 

which is updated by measuring the incident numbers I\ t2 and employing Eq. ([51]). 

With the distribution at hand, we could proceed similarly as in the previous Section 3.1 , try- 
ing to estimate the average path length related to the increasing incident numbers in particular, 
in order to gain a measure of the change taking place in the system. 

However, in the following we proceed differently, in a way which appears more suitable 
to further generalization. We introduce the Fisher-Rao information metric for the purpose of 
quantifying the change due to the chaotic, even if deterministic motion from one configuration 
to the next [[0], 12]. It is the uniquely determined Riemannian metric (except for an overall 
multiplicative constant) on the space of states which are probability distributions [fTOfl . In our 
present case the states are simply described by the pair of coordinates Q l = (Xj) \i=i, 2 G [0, 1], 
considering the coordinate H int to be fixed at a constant value. Then, the 'distance' ds between 
the states P(x\H int ; Q k ) and P(x\H int ; Q k + dQ k ) is given by: 



ds 2 = gijdQ'dQi 



with the metric: 

9ij 



J dxP(x\H mt ;Q k 



d\nP(x\H int] Q k ) d\nP{x\H int] Q h 



J(^,e)- J(e,e) 

J(M) 



(Q 



(59) 



(60) 



(61) 



employing Eqs. (||0|), ([42|), ([43[), and particularly Eq. ( p8[) in the explicit calculation. Note that 
the coordinates, being probabilities, are constrained by Q 1 + Q 2 = 1, which further simplifies 
the result obtained here. 
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However, we would like to express Eqs. (p9|)-(|61|) in terms of the directly measured incident 
numbers. Thus, with the help of Eqs. fl5T|), we obtain: 



As = 2 J^)-J^e)\I 2 ^- h AI 2 \ 



J(l,l) y/hh(h + L 



2 



^ fe.t.y , (63) 



J(l, l)(/i + I 2 ) \V A \M2 

where the second equality is due to the dichotomous character of the pair of increments AIi and 
AI 2 ; i-e., we have AI^ = 0, 1, however, the two increments present the results of two mutually 
exclusive measurements which cannot happen in coincidence. 

Several remarks are in order here. First, the Eq. ( |63"D is invariant under the similarity 
transformation I*. — > £ • Ik \k=i,2 and As — > As/£. Thus, generally, the amount of "change per 
incident" decreases inversely with the total number of incidents, i.e., loosely speaking, with the 
age of the evolution going on. This seems natural in a statistical context, but less familiar in 
the context of dynamics, where we are accustomed to a linearly flowing time. 

Second, the dichotomous increments of the incident numbers can be mapped on a binary 
sequence of zeros and ones, for example, according to (Aii = 1) — >• and (A/2 = 1) — ► 1. Then, 
the number of zeros and ones in the produced bitstring is simply related to I\ — 1 and I 2 — 1, 
respectively. Thus, one could compare the behavior of the function As with corresponding ones 
for different dynamical models or with randomly generated sequences. 

Third, the statistical factors in the results of Eqs. (|6T|) (|62|) , encoding the dynamics, are 
nicely factorized from the quantities incorporating the evolving observables I\ and I 2 ; an anal- 
ogous factorization is seen in Eq. fl35l). We believe that the very simple structure of our results 
is indeed of a more general kind and similarly will be encountered, whenever one succeeds to 
identify pairs of dichotomous observables, which somehow reflect the symmetry of the system. 

Furthermore, the generalization involving n-tuples of exclusive observables can be worked 
out along these lines. 

Coming back to the construction of a reparametrization invariant "time" T, we now set: 

AT.^l, (64) 

with N = ii + I 2 , i.e. stretching the time with an extra factor N, in order to compensate for 
the aging of the measure As, which we mentioned. Furthermore, we divide out the square- 
root-of-energy factor, which is determined by the initial conditions. This results in a scaling 
with energy, similarly as in Section 3.1. At high energy H int , Eq. equals (two times) 

the kinetic energy, which might seem more appropriate here. However, it is by itself not a 
reparametrization invariant quantity in the present interacting case. 

We calculate the "time" passed by accumulating the time steps AT. The constant, though 
energy dependent, overall prefactor contained in the result of Eq. (^4j) is set equal to one, re- 
calling that the Fisher- Rao metric which enters here, see Eqs. (|59|)-(|6"T|), is unique only up to 
a constant factor. The Fig. 4 shows the resulting behavior of the parametrization invariant 



"time" T as a function of the proper time r, Eq. (PB[), which is also incorporated in the present 
numerical simulations. For the two initial conditions, differing in the initial velocities, we ob- 
serve the expected scaling with energy. Furthermore, while the size of individual fluctuations is 
comparable to what was found in the noninteracting case of Fig. 2, we notice here characteristic 
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deviations from the average linear dependence over longer proper time scales. Such nonlinear 
behavior becomes even more prominent in Fig. 5, with long time excursions. 

In any case, this demonstrates the existence of a monotonous "time" function also for the 
present more strongly irregular system, as compared to the ergodic (but integrable) one of 
Section 3.1. Presumably, an approximately linear dependence will only become recognizable 
for much longer runs, i.e. for larger incident numbers. On the other hand, of course, the 
statistical construction of T based on quasi-local observables could be refined by considering a 
higher-dimensional state space from the outset, cf. Eq. fl3~7D , i.e. incorporating more observables. 

This completes our discussion of two examples of the construction of a coarse-grained 
parametrization invariant evolution parameter. 



4 Stroboscopic quantization? 

It seems worth while to draw attention to a possible connection between the present subject 
and the deterministic classical systems which are effectively quantum mechanical, as recently 



discussed by 'tHooft, see Ref. |13| and further references therein. 

The total number of incidents obtained from our exclusive observables always increases by 
one. These are the "ticks of the clock", numbered by n e N. However, via the constructed 
"time" sequence {T n }, the proper times r(T n ) vary irregularly, in general, due to the quasi- 
periodic internal motion. Therefore, the free external motion, Q\\(T) = q\\(r(T)) = +r(T)p, 
shows stochastic behavior in terms of T. We emphasize that r(T) here is not directly related 
to r(t) introduced in Eq. (|36|). 

Let us assume that we could exhaust all possible invariant quasi-local observables and cannot 
further resolve the "time between the incidents" marking T and T + ST. 

That there arises a discreteness in the optimal "time" sequence generated in this way in 
a generic nonintegrable system |J seems plausible for two reasons. First, by realistically re- 
stricting the observables to be localized near to the observer, we loose information about the 
full trajectories in coordinate space. Second, by insisting on gauge invariant observables, we 
loose the possibility to make use of additional momentum space observables (except for a pos- 
sible limited number of constants of motion), or higher proper time derivatives, in order to 
reconstruct the full trajectories and a smooth evolution parameter based on them. 

Anyhow, for given initial conditions, the incident counts as well as the coordinates q\\(T n ) 
present discrete physical data which are deterministically determined in our model. (The 
constant orthogonal components are of no interest here.) Then, the question arises, can we 
still give an invariant meaning to and predict the future of the free external motion, instead 
of letting the system evolve and measuring? - In general, the answer is 'No', since neither the 
future discrete values of T, see Eq. fl64|), for example, nor the value of r at those future instants 
will be predictable. 

However, motivated by the two practical examples of Section 3 , we introduce the notion of 
a well-behaved "time", which will be sufficient to maintain some predictability indeed. 

Let a well-behaved "time" T be such that the onedimensional motion of a free particle, co- 
ordinated in a reparametrization-invariant way by the sequence {q(T n )}, is of limited variation: 

{T n } is well-behaved 

<^=> 3 f , At : fn - At < r(T„) < fn + At , Vn G N 
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3 q, Aq : qn - Aq + g (0) < q(T n ) < qn + Aq + <? (0) , Vn G N 



(65) 



with g = fp(°), Aq = Arp(°\ and where q(°> ,p(°> denote the initial data. This is the next-best 
we can expect from a reasonable "time" , motivated by Newtonian physics here. It is important 
to realize that there is no hope for continuity or periodicity, since our construction of T is based 
generally on the ergodicity of a quasi-periodic process and related Poincare sections. 

Assuming a well-behaved "time", we see that the sequence of points {q(T n )} can be mapped 
into a regular lattice of possibly overlapping cells of size 2 Aq and spacing q. "As the clock ticks" , 
i.e. n — > n + 1, the particle moves from one cell to the next. On the average this takes a proper 
time f and physical "time" T. Here we assume for convenience that {T n } itself is of limited 
variation, 3 T, AT : T n = Tn ± AT, Vn e N. In the simplest example, in Section 3.1, we have 
f = S~ 1 (7Ti +7r2°^) _1 and AT = 0, showing a dependence on the initial momenta and the scale 
of localization of the observed incidents. 

Thus, for a well-behaved "time" sequence (of limited variation), we have mappings from 
"clock ticks" to ( "time" as well as) space intervals. We identify the space intervals, containing 
the respective q(T n ), with primordial states: 



\n 



= [qn — Aq + q^°\ qn + Aq + q 



,(o)l 



(66) 



In order to handle the subtle limit of an infinite system, we impose reflecting boundary condi- 
tions and distinguish left- and right-moving states. Finally, then, the evolution is described by 
the deterministic rules: 



n 



\n 



| AO 



n+ 1 , neN , 

\n + 1) , -iV + 1 < n < N - 1 

l-N + l) , 



(67) 
(68) 
(69) 



with 2N states in all, | — N + 1), ... , \N); states with a negative (positive) label correspond 
to left- (right-) moving states, according to negative (positive) particle momentum. Reaching 
the states |0) or \N), the particle changes its direction of motion. 

The rules (|67D-(|S9|) can altogether be represented by a unitary evolution operator: 



U(5T = T) = e 



-iHT 



-in(N+l/2)/N 



( o 

1 

1 
V 



l \ 







(70) 



which acts on the 2iV-dimensional vector composed of the primordial states; the overall phase 
factor is introduced for later convenience. Here T is the natural scale for the Hamiltonian H. 

In the following analysis we apply 't Hooft's method, who considered the discrete and strictly 
periodic motion of a classical particle on a circle, introducing a Hilbert space based on the 
primordial states fl3|]. The evolution operator turns out to be diagonal with respect to the 



discrete Fourier transforms of the states \n). We define the basis functions: 



{k\n) = f k (ri) 



(2iV)- 1/2 exp 



iitkn^ 

■ N ' 



(71) 



Time without time: a stochastic clock model 



15 



They present a complete orthonormal basis, since Y<k=-N+i( n '\k) (k\ n ) — &n>ni with [n\k) = 
(k\n)*, and noting that fk{n) = f n (k). This yields indeed: 

1 2%i 1 

(k'\U\k) =5 k/k exp(-m[N+--k}/N) & (m'\U\m) = 5 m , m ex P (-^—^[s + m+ -]) , (72) 

where the equivalence follows from 2s + 1 = 2N and relabeling and replacing the states \k) by 
states |m), with m = —k + 1/2 and — s < m < s. The phase factor of Eq. ( |T0[ ) contributes the 
additional terms which assure a positive definite (bounded) spectrum similar to the harmonic 
oscillator case. 

Recalling the algebra of SU(2) generators, and with S z \m) = m\m) in particular, we obtain 
the Hamiltonian: 

H = , 2n ,- (S z + s + -) , (73) 
(2s + 1)T V 2 J ' v ; 

i.e., diagonal with respect to |m)-states of the half-integer representations determined by s. 

Continuing with standard notation, we have S 2 = S 2 + S 2 + S 2 = s(s + l), which suffices 
to obtain the following identity: 

Furthermore, using S± = S x ±iS y , we introduce coordinate and conjugate momentum operators: 

q=\(aS^ + a*S + ) , p = l -(bS„ + b*S + ) , (75) 

where a and b are complex coefficients. Calculating the basic commutator with the help of 
[S + , S-] = 2S Z and using Eq. ([73]), we obtain: 

[q,p] = i{l--H) , (76) 

7T 



S 2 + S 2 = (2S j 1)2 (|o| V + \b\ 2 q 2 - (3a -Sb + Vta- %6){g,p}) . (77) 



provided we choose ^(a*^) = — 2(2s + 1) . Incorporating this choice, we calculate: 

4 

In order to obtain a reasonably simple Hamiltonian in the infinite system limit (s — > oo), we 
finally choose: 

ff , 2 Pk 

(78) 



7T ' 2s + 1 V T 

Then, defining = 2n/(2s + 1)T, the previous Eq. ( j74|) becomes: 

= lp2 + 1^2 2 + £ ( 1 2 + ^2) (7Q) 

2 F 2 y 2tt 4 ; ' v ; 

showing a nonlinearly modified harmonic oscillator Hamiltonian at this stage, which corre- 
sponds to the result of Ref . [|TB| . 

Following our construction of the reparametrization-invariant "time", we presently keep T 
finite while considering the infinite system limit. Thus, for s — > oo, we have uo — > and obtain: 
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which has the low-energy limit H m p 2 /2. On the other side, the energy is bounded from above 
by ir/T, since we have (T/ti)p 2 = 4(2s + l) -2 ^ 2 (< 1, for s — > oo, when diagonalized) . 

Interestingly, towards the upper bound the violation of the basic quantum mechanical com- 
mutation relation, as calculated in Eq. ([76]), becomes maximal, i.e., there q and p commute and 
behave like classical observables. 

We conclude here with calculating the matrix elements of the operators q and p between 
primordial states. Relating the primordial states \n) to the |m)-states with the help of Eq. (|71~D 
and using S±\m) = [s(s + 1) — m(m ± 1)] ly/2 |rri ± 1), we obtain: 



2s + 

s 

^ exp(g T (n' -n)(m-l/2))[s(s + l)-m(m + l)] 1 / 2 . (81) 



For large s the summation is replaced by an integral. Then, using Eqs. ( |75|) and (|78|), the results 
are (s — > oo): 

(n'Mn) = l^ M<n'-n))n' + n 

y m ! 2 n'-n 2 K J 

. . . . . 1 IT Ji(ir(n' — n)) , . 

(n'pn) = ~ , 83 

2 V vr n — n 

where J\ denotes an ordinary Bessel function of the first kind. Thus, neither the position nor 
the momentum operator is diagonal in the basis of the primordial states. 

In any case, we find here once more an emergent quantum model based on a deterministic 
classical evolution, similar as in Ref. fL3 |. 



The interesting feature presently is that the construction of a reparametrization-invariant 
"time" based on quasi-local observables naturally induces the particular stochastic features in 
the behavior of the external particle motion. The remaining predictable aspects of its motion 
can then most simply be described by a unitary discrete-time quantum mechanical evolution. 



5 Conclusions 

Our construction of a reparametrization-invariant "time" is motivated by the observation that 
"time passes" when there is an observable change, which is localized with the observer. More 
precisely, necessary are incidents, i.e. observable unit changes, which are recorded, and from 
which invariant quantities characterizing the change of the evolving system can be derived. 

A basic ingredient is the assumption of ergodicity, such that the system explores dynamically 
the whole allowed energy surface in phase space. Generally, then there will be sufficiently 
frequent usable incidents next to the observer. 

We illustrate this in Section 3 with a simple timeless model of a nonrelativistic particle 
moving classically in five dimensions, of which two internal ones are compactified to form a 
torus. Defining suitable window functions, localized observations of incidents correspond here 
to the trajectory passing by the window. Thus, the incidents reflect properties of the dynamics 
with respect to (subsets of) Poincare sections. Roughly, the passing "time" corresponds to the 
observable change there. 
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In the first example in Section 3.1 the dynamics is very simple and even integrable ||. We 
know that the system is ergodic in this case and use this fact, in order to calculate statistically, 
based on quasi-local incident counts, the invariant path length run by the system. This provides 
an invariant measure of "time" . We find that the proper time is linearly related to this, however, 
subject to stochastic fluctuations. 

In the second example, Section 3.2, we have additional interactions which mix the trajecto- 
ries in a much more complicated way. We assume that the system is ergodic. Our statistical 
formalism then allows to measure the change of the phase space distribution which is updated 
by the information coming from the incident counts. Essential ingredient is the Fisher-Rao 



information metric [10, 11, 12 1. This calculation again provides an invariant "time" function, 
which behaves qualitatively similar as in the previous example, however, with characteristic 
long-time deviations from an average linear relation to proper time. 

It will be interesting to generalize the present model to relativistic systems, as well as to 
generalize the statistical formalism. A related question is, which are the (discrete) limitations 
of such a construction of "time" based on localized observables, i.e. how close can one get to a 
linearly flowing Newtonian time, in a generic nonintegrable system. 

Finally, we pointed out in Section 4 a possible connection of the present subject to the study 
of quantum mechanical systems which have an underlying deterministic classical dynamical 



model ||13|| . Indeed, we find that for certain "well-behaved time sequences" the remaining 
deterministic aspects of the induced stochastic classical motion of the (in our model) external 
particle can be most simply described by a discrete-time quantum mechanical Schrodinger 
evolution. This "stroboscopic quantization" may arise under more general circumstances, if 
one insists to construct by statistical means an invariant "time" from localized observables. 
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Figures: 

Fig. 1: A typical trajectory. Initial conditions for the free motion on the torus are <j>^\ = and 
7r[°' ) = 1, 71*2 = y/3 (arbitrary units). Here t — r (proper time); the final time is r = 50. 

Fig. 2: The reparametrization invariant "time" T, Eq. fl35|), based on localized incident counting, 
as a function of (computer) proper time r for two different proper time intervals, cases a) 
and b); initial conditions as in Fig. 1; window parameter: 5 = 1/3. The lines result from 
linear fits. 

Fig. 3: A typical trajectory with interaction parameters u — Q — 3. The initial conditions are: 
0i°2 = an d 7r|°' ) = 3, 7i"2 — \/5; the final proper time (t = r) is r — 30. 

Fig. 4: The reparametrization invariant "time" T, from Eq. ( |64"D with constant prefactor set to 
one, as a function of (computer) proper time r for two different sets of initial conditions: 
^i°2 = (always), = 5, 712 = \/l7 (upper curve) and tt^ = 3, 712 = v^5 (lower 
curve); interaction parameters: u — — 3, window parameter: 5 = 1/3. Lines from 
linear fits. 

Fig. 5: Same as Fig. 4 (upper curve), illustrating nonlinear long-time behavior. 
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